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We use the coherent state approximation for the skyrmion in 
the isotropic quantum Heisenberg ferromagnet to obtain analytical 
expression for skyrmion hopping integral as a function of distance. 
It appears the skyrmion hopping is restricted by distances larger 
than its effective diameter. Some puzzling properties of the orbital 
skyrmionic modes are discussed. Spin distribution in these states 
is obtained and effect of quantum spin contraction is demonstrated. 
The principal possibility of lowering the skyrmion energy due to a 
bonding-like state formation is illustrated. 
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1 Introduction 

Translational motion of the 2D topological excitations, considered 
as a particle-like object, is a problem of enhanced importance both 
in theoretical and practical sense. Classical approach implies mak- 
ing use of phenomenological equations of motion. A variational ap- 
proach is one of the most attractive methods to study an appropri- 
ate quantum problem. The variational wave function of collective 
motions, like translations, rotations and oscillations, was initially 
proposed by Hill and Wheeler |J to be 

*([r]) = / rf 2 R/(R)vM[r]), (1) 

where \&R,([r]) is a particle-like wave function for the topological 
excitation centered at point R = (R, 9), /(R) a weighting function, 
and [r] is a set of particle coordinates. For continuous 2D media the 
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angular dependence of the weighting function follows immediately 
from the cylindrical symmetry of the system, so, 

f(R) = f l (R)exp(iie) 

with / being the angular momentum. By varying the weighting 
function to minimize the expectation value of the Hamiltonian 

£=(*([r])|#|*([r])>, 
one obtains an integral equation for /(R), 

J d 2 K'K{K', R)/(R') = 0, (2) 

where the kernel is 

K(R', R) = (* w {[r])\H - E\* n ([r})) = H(R', R) - ES(R', R), 

(3) 

where we have made use of matrix elements of Hamiltonian, and 
overlap integral S^R', R). The kernel K(R',R) may be considered 
as an effective transfer (hopping) integral for quasiparticle transfer 
between R' and R sites. 

It should be noted however that full realization of above men- 
tioned scheme (|]||) for the real topological defect on the lattice is 
a difficult task as the lattice has no authentic rotational symmetry 
unlikely the conventional case of the vortex in the superfluid [0. Ob- 
viously, the first step to treat this problem implies using continual 
approximation for describing soliton in the lattice. 

Below we calculate quantity K(R', R) for quantized skyrmion 
being peculiar topological defect described in the framework of the 
coherent state approximation and use continual approximation to 
treat skyrmion motion in a discrete lattice. 

Skyrmions are general static solutions of 2D Heisenberg ferro- 
magnet with isotropic spin-Hamiltonian 

H = — JSj^SjSj+b, (4) 

obtained by Belavin and Polyakov || from classical nonlinear sigma 
model. A renewed interest to these unconventional spin textures 
is stimulated by high-T c problem in doped quasi-2D-cuprates and 
quantum Hall effect. 
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The skyrmion spin texture consists of a vortex-like arrangement 
of the in-plane components of spin with the ^-component reversed 
in the centre of the skyrmion and gradually increasing to match the 
homogeneous background at infinity. The spin distribution within 
classical skyrmion is given as follows 

2rA 2rA . r 2 - A 2 

S - = ^^ cos ^ S v = -2-^ sm( P> Sz = ^^- ^ 

In terms of the stereographic variables the skyrmion with radius 
A and phase (po centered at a point zq is identified with spin dis- 
tribution w(z) = jz^, where z = x + iy = re lip is a point in 

the complex plane, A = \e t6 , and characterized by three modes: 
translational, or positional ZQ-mode, "rotational" ^-mode and "di- 
latational" A-mode. Each of them corresponds to certain symmetry 
of the classical skyrmion configuration. For example, #-mode cor- 
responds to combination of rotational symmetry and internal U(l) 
transformation. 

The simplest wave function of the quantum spin system, which 
corresponds to classical skyrmion, is a product of spin coherent 
states H]. In case of spin s = ~ 

<M0) = E[cos 6 ^ | T ) + sin|e-^ |j>], (6) 

r 2 -A 2 

where 6i = arccos r \ +X 2 ■ Coherent state implies a maximal equiva- 
lence to classical state with minimal uncertainty of spin components. 

Classical skyrmions with different phases and radia have equal 
energy. Nevertheless, stationary state of quantum skyrmion is not a 
superposition of states with different phase and radius || and has a 
certain distinct value of A. Fluctuations in the in-plane orientation 
and size A for quantum skyrmion are intimately connected; rotating 
a skyrmion changes its size. This result follows from the commuta- 
tion relations of quantum angular momentum operators 0. 

In our recent paper it has been obtained an analytical expres- 
sion for overlap integral of skyrmion states 



S(Ri, R 2 ) = exp[-^-(i? 2 2 - 4zA 12 )], (7) 
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if R 12 < 2A or 



5(R l5 R 2 ) = exp[-^-(i? 2 2 - AtA 
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,i?i2 + A/i??2-4A 2 



- R^Rh - 4A 2 + 2A 2 l < 12 V 12 = )], (8) 

Hia - xAr?, - 4A 2 



Li2 - y-n-12 

if R±2 > 2A. Here one introduces a quantity 

A12 = ^ I -Ri || -R 2 I Sin((y9i - </? 2 ) 

being area of a sector skyrmion covers while moving in the plane, a 
is a lattice parameter. Skyrmion, having wandered closed contour 
on the plane, acquires the phase Att^ 1 , that is each spin lying inside 
the contour contributes An into the total phase of the skyrmion, so 
its motion appears to be likely the motion of a charged particle in 
magnetic field || || [UJ . Extending this analogy one may introduce 
the effective magnetic length d = (^r)~ 5 . Note the dependence of 
the phase factor in the overlap integral on the origin of the vectors 
Ri, R 2 . It should be noted that for large separation i? i2 2> A the 
R dependence of the overlap integral obeys the power law 

5(R 1 ,R 2 )oc(A/ J R 12 r 

with p = (X/d) 2 . 



2 Hopping integral for stiff quantum skyrmion 

One of the main goals of the paper is a calculation of the nondiagonal 
matrix element of the Hamiltonian (Q) 

E{R 1) R a ) = «([r])|E|<([r])> ) 

between the states of quantum skyrmions, located at points Ri and 
i?2, respectively, or transfer integral, likely being the principal quan- 
tity defining the quasiparticle motion of this topological defect. 

We use approximation of "stiff" skyrmions, which takes no ac- 
count of internal skyrmion quantum fluctuations and is closely re- 
lated to classical limit. Making use of the coherent states (Bl) we 
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write out the appropriate matrix element as follows: 







Q+fc) (<M+i?12 













i/(Ri,R 2 ) — — -JSizYiij 
For the spin coherent states the following relations exist : 



(q\S z \fi) = -S- 



(q I S+ I /i) = 2S 1 - 



(f 1 


5+ 


/i) 




/i) 



1+U| 2 ) S (1+|//| 2 ) S ' 

i+M 2 )s(i+|/i| 2 )s' 

s 1 -^, do) 

1 + q^ 

2S— ^— , ^1 V ' r/ = 2,S'- 





5, 




( 




//) 





s 






(f 1 


/i) 



?/i I /i) 1 + 

Taking account of the spin distribution in the skyrmion we write 
out the matrix element addressed as follows 

H(R U R 2 ) = ~ JS 2 S 12 £^[[((q - - Rz) ~ A 2 )(fe +b - x 

x (q^~ b - R~ 2 ) - A 2 ) + 2 A 2 (q +6 - fli ) (q - R~ 2 ) + 2 A 2 - TQ (s - 
[(A 2 + (s - i20(ft - ^))(A 2 + - - R2))]- 1 ]. (11) 

Below we are going to find the limit of the right hand side when 
a — ► with the constant a of lattice assumed to be square. To do 
it one has to make expansion in lattice constant a and to keep the 
three first terms, then to turn the sum into integral. The first (zero 
order in a) term in the expansion, as it is easy to see, is simply the 
ground state energy. The linear in a term in the expansion is zero 
due to the inversion symmetry. So, one has to define the second 
order term for a given site i. Straightforward calculation gives for 
that second order term in a in the matrix element for scalar product 
S fo SjS i+ fe the following result 

-8A 2 2 
a. 



(\* + (q - R 2 )(q - Rj)* 
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While coming to the limit a — > the summation in ( |TT| ) is replaced 
by the integration, and this term results in a nonvanishing contri- 
bution to the transfer integral: 

/•CO 

HCRx, R 2 ) = [E + US 2 \ 2 / rdr 

Jo 

(A^ + (re lv — R 2 ){re 1 ^ — Ri)) 

where En is a ground state energy of spin system. An angular inte- 
gral is given by following expression 

dip- 



I 



o (A 2 + {re~ iip - R 2 )(re i v - Ri)) 2 

zdz 2n(xi + x 2 ) 



%{z — Xi) 2 (z — x 2 ) 2 (rR 2 ) 2 (rR 2 ) 2 (xi — x 2 ) 1 
so the radial integral is 

2nr(r 2 + Rji 2 + X 2 )dr 1 



r 

Jo 



\(r 2 + RiR 2 + A 2 ) 2 - 4,r 2 R x R 2 f A ^ 

where 

x (rR R) * + r 2 + R 1 R 2 ±K 

Xi, 2 {r,Ki,K 2 ) — _ . 

2rRo 



K{r, R Xl R 2 ) = y/(A 2 + r 2 + R 1 R 2 ) 2 - Ar 2 R ± R 2 , 

where the complex root has positive real part. 

When i?i2 < 2A we finally obtain for the transfer integral 

#(Ri, R 2 ) = [E + 4tt J5 2 ]5(Ri, R 2 ). (13) 

As R 12 = ff(Rx,Ri) = [E + AttJS 2 }, and this is the classical 
skyrmion energy. The structure of the transfer integral if(Ri,R 2 ) 
is such, that its contribution to the effective hopping integral (3) 
K(Hi, R 2 ) at R\ 2 < 2 A is exactly compensated by the overlap con- 
tribution. Thus, effective hopping integral for the "stiff" skyrmion 
is zero for distances R\ 2 < 2A. 
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When R12 > 2A, the integral (|12|) contains two singular points 
and, respectively, it is ill-defined in that region. To avoid the con- 
tribution from these points and obtain result independent of choice 
of coordinates we calculate its principal value, that is remove the 
contribution of points, that lie inside the circles of small radius e 
with centers in singular points. To this end we divide the plane into 
two equal parts drawing the border between halfplanes at the mid- 
dle of the line connecting two singular points. The value of integral 
( fi"2|) is twice the value of integral over halfplane. Integral (|l"2"D being 
represented as halfplane integral in coordinate system with center 
in singular point takes the form 

roc 

H{R h R 2 ) = [E + 8JS 2 X 2 J rdrx 

/•2tt— arccos — j 

X Arccos s d(f> r 2 (r + - c)e~^ - (4^ + c)e^) 2 + ' ' ' ' 

+8JS 2 X 2 l\dr dtp— -5 — ^ — — }S(R 1 ,R 2 ), 

Jo Jo Y r 2 (r + -c)e-^ - + c)e^) 2J V ' ; ' 



where c = y ^p- — A 2 is a half-separation between singular points. 
Here we have restricted ourselves to the case of real Ri, R 2 to sim- 
plify calculations. It is easy to see that last term in (O) is zero as 
angular integral is zero no matter what the value of r is. Then, 

d<f (r + (4^ - c)e-*v ~{ B f + c)e^) 2 ~ ^ 

1 r e** e iv a;_ + x +1 e^-ac., 

m ■ 



%{r 2 + AX 2 ) e l(p — x_ e lLp — x + X- — x + e %ip — x + 
where 



x +, 



± Vr 2 + AX 2 



2{^f + c) 

As a result, we obtain the following radial integral 

roo r2it— arccos - | 

/ rdr / r dip- 



arccos 2 r 2 (r + - c)e-**» - + c)e^) 2 
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%r 2 + 4A 2 ) [ J (r 2 + AX 2 ) X (16) 



x arctan \ (r 2 + 4A 2 ) (— -) — — - / = 0. 

V c 2 r 2 -4c 2 A 2 + r 2 (r 2 + 4A 2 ) 

Thus, it appears that the effective hopping integral is nonzero 
in this regime and its value is simply proportional to the overlap 
integral. 

So, the final result for effective hopping integral is 

K(Rx,R 2 ) = -47r5 2 J5(Ri,R 2 )^(i? 12 - 2A), (17) 

where 

0(i? 12 -2A) = (!' ' li i ? 12>2X , 
v 12 ; \0, otherwise ' 

and 5 , (Ri,R 2 ) is given by expressions (|7|), (§). 

Thus, one might conclude that the skyrmion hopping is restricted 
by distances larger than 2A. The nonzero effective hopping integral 
for distances larger than 2A allows to lower skyrmion energy by 
construction of appropriate combination of skyrmion states centered 
at different points of lattice and having definite value of angular 
moment. Energy of these states will depend on the value of the 
angular moment and on number of the levels as it is the case for the 
Landau levels in a magnetic field. Note also that transfer integral 
has the negative sign as must be the case. 

The states considered have definite phase and radius. In fact 
phase and radius can not be treated separately as radius is deter- 
mined by spin of the skyrmion and spin in turn is canonically conju- 
gate to the phase. It is easy to show that states of definite phase and 
radius are not mixed by the Hamiltonian. This is related to the fact 
that we deal with the skyrmion of charge q = 1, which has infinite 
spin. As the spin and phase operators do not commute, the spin 
must have infinite value for the phase to be definite. In the case of 
higher topological charges there would be mixing of the states with 
different phases and radius. The similar conclusions were drawn by 
Stern 0. 
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3 Orbital motion of skyrmions. Small and large 
orbital /-skyrmions 

Making use of the above derived results we may exploit particu- 
lar analogy between skyrmion and particle in zero Landau level to 
construct a subset of orthogonal states with definite value of or- 
bital momentum, or /-skyrmions. Consider firstly the following or- 
bital skyrmionic modes being linear combinations of single skyrmion 
states 

4f (R) = J\xv (iimt^dQ, (18) 

where ^^ exp ( i0 ) * s a single skyrmion state centered at the point Z = 
Rexp(iQ). Below, we will distinguish small and large /-skyrmions 
depending on whether R < A or R > A, respectively. 



3.1 Small orbital /-skyrmions 

It seems, any two equally centered ^f k (Ri) and ^f k (R 2 ) states with 
commonly different R\ < A and R2 < A, but the same orbital mo- 
mentum / are linearly dependent. Indeed, for small /-skyrmions one 
might use a simple Gaussian form (|7|) for the overlap integral of 
single skyrmionic states, so the corresponding unnormalized overlap 
integral is given by a rather simple analytical expression 

S(R X , R 2 ) = / exp[i/(9 1 - a )] exp[- — {R\ + R\ - 2R,R 2 
Jo 4cr 

«p[-i(ei - e 2 )])]de x de 2 = ^ exp[-i(j2? + Rl)}{^) 1 . 

(19) 

It is easy to see that after normalization the overlap integral for two 
/-skyrmions with different R\ and R 2 

S(Ri,R 2 ) 



S(Ri, Ri)S(R 2 , R 2 ) 



does not depend on Ri, R 2 at all, which means at least that the 
corresponding functions are not independent. So, these states are 
distinguished by a single quantum number / which is an orbital 



momentum and are degenerate in energy in full analogy with a zero 
Landau level. It should be emphasized that for "/-skyrmion" \E f f fc (i?) 
with R < A the orbital moment is restricted by non-negative values 
I > 0, due to the definite sign of the Berry phase in the overlap 
integral for single skyrmion. To describe a spatial distribution of a 
" single skyrmionic (quasiparticle) density" for the orbital /-skyrmion 
one might use an overlap integral of the orbital state (|18D with 
the simple skyrmion state in the point z = rexpiO which has a 
rather simple form for r < A 

S l sk (r,R) = J\xp[iie] exp[-^(r 2 +R 2 -2rRexp[-t(9-e)])]de = 

2tt 1 r /? 

- exp[z/0] exp[-^(r 2 + R 2 )}{^)\ (20) 

that with account for normalization conditions gives the quantity 

sL(rH ^||^ exp|!W,( ^ , ' exp[ -^ ] ' (21) 



which does not depend at all on R value. This distribution function 
has maximum at r = r\, where 

so, the orbital momentum I defines the characteristic "orbital area" 
nrf measured in terms of a "flux quantum" This relation looks 
like a specific quantization condition for the "orbital area". The T\ 
value defined by expression (|22"D may be called as a characteristic 
"internal" radius of the orbital states which magnitude does 

not depend on R value. One should note that for "s-skyrmion" 
r = as it should be for typical quasiparticle orbital s-states. 

Actually, the maximal value of orbital momentum for the small 
/-skyrmion is restricted by the size of the bare skyrmion and, in 
general, can not significantly exceed l\, because of the number of 
translational degrees of freedom is restricted by the "skyrmionic 
area" ttX 2 . If we consider values of / much larger than l\ we would 
obtain meaningless result for the spin distribution in the /-skyrmion 
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because in this case continual approximation for the skyrmion over- 
lap breakes down and expression (|7|) that neglects discretness effects 
is not adequate. 



3.2 Quantum spin contraction for small radius /-skyrmions 

Spin distribution in the puzzling "Z-skyrmionic" ^i(R) states with 
small R < A can also be shown to be independent on radius R. Be- 
sides, the corresponding spin texture is characterized by a quantum 
spin contraction, that is the mean spin length is less than S. It can 
be calculated with making use of relations (|TUD, so that we have for 
the spin averages 

1 f 27T I (*~*0 I c I (£z£2)\ 

s+(z) = -J Q ex P [uie, - e 2 p ^ ( J zi) ±^- ^l Sl2 ( Zl -z 2 )de 1 de, 



A I A 



1 f»* vua amc 2A(z-rexp(i0 1 )) 

= — / exp \ll{0x - 2 ) St7 7 — 7T7TTV7 ~ ; rrrX 

N Jo L V n \ 2 + {z-rexp(i6 1 ))(z-rexp{-i6 2 )) 

x exp[-^(r 2 - r 2 exp \-i(Q x - 9 2 )])]dB 1 d9 2> (23) 

where z\ = r exp z 2 = r exp i9 2 are locations of distinct skyrmion 
states on the circle of radius r, 

(2it) 2 r 2 r 2 
N = S(r,r)= { —L(—YeM 



l\ y 2d 2J 2d 21 

is a normalization integral. Introducing complex variables i\ = 
exp(i^i), z 2 = exp(— z# 2 ) and taking integral in £\, and then in z 2 we 
obtain 

S ^) = ^l I X2 f KZ ~ , , ^{-^y-rH^dzxdz, 

N J J A z + [z — rz\)[z — rz 2 ) 2d A 

(24) 

l,2d\, , d' d l . 2X(z-rz l ) . r 2 , , .. 

= Try— 7t) hm ——r hm — r — ; — — —exp — ~Z\Z 2 \\. 

l\ K r 2 ' z 2 -,o dz l 2 £i-*o dz\ L A 2 + (z - rz x ){z - rz 2 ) l 2d 2 JJ 



This is an exact expression as integrals in (|24"D have no other poles 
but those of (I + l)-order in i\ = and z 2 = 0, as long as r < I. For 
the S z we have a similar expression 
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l\ r 2 z 2 ^o dz\ *i-»o dz\ X 2 + (z — rz x ){z — rz 2 ) 2d 2 

(25) 

Surprisingly, for the case I = we obtain the same spin distribution 
as in the simple single skyrmion state. In other words, in the frame- 
work of a classical approach the s-skyrmion is equivalent to a single 
skyrmion. For the "p-skyrmion" with I = 1 the straightforward 
calculation gives 

2Xz 8Xd 2 z 8Xd 2 z 2 z 2Xz . 4d 2 X 2 , 



A 2 + 2z (A 2 + ^) 2 (A 2 + A 2 + 22 (A 2 + zz) 21 

z = T^— - 1 " 7T2— 26 

A 2 + zz (A 2 + zz) 1 

where all values are measured in lattice constant. Here, we deal 
with a spin length contraction, compared with the case of classical 
skyrmion, specified by an universal function of skyrmion radius and 

r =| z | 

4d 2 A 2 , 



F(A,r) = [l 



(X 2 + zz) 2> 

Increasing the skyrmion radius entails decreasing the spin contrac- 
tion. Minimal spin contraction takes place for large r's and A's, 
while the maximal spin contraction does in the center of the spin 
texture for small A's. Interestingly, that for small skyrmion with 
X = 2d the effective spin length in the core of the orbital p-skyrmion 
turns to zero: S(r = 0) = 0. However, as it is easy to see for very 



small bare skyrmions with A < c( the expression (26) leads to un- 
physical results. Indeed, the continuous approximation used above 
appears to be irrelevant for small skyrmions with radius A less than 
the effective magnetic length d. 

The case / > 1 can not be represented in the simple form like 
(|26|), however, with increase of I the value of spin contraction in the 



center of spin distribution also increases. It should be noted, that for 
the orbital skyrmions with I — l\— 2n \ = (which is just a ratio 
of "skyrmionic area" 7rA 2 to the area corresponding to the quantum 
of flux t^), the mean spin length in the center of spin distribution 
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turns to zero, that is the spin contraction appears to be maximal. 
This is due to the fact that at I = l\ the characteristic internal 
radius r; = / A of the orbital skyrmion defined by expression (|2"2"|) 
coincides with skyrmion radius A, so in the /^-skyrmion core we 
deal with a superposition of vectors with zero z-projection and all 
possible phases. 



3.3 Large radius orbital skyrmions 

Above we addressed to the particular type of the small orbital 
/-skyrmion with bare radius R < A. Unfortunately, for large / 
strongly exceeding l\ continual approximation breakes down so we 
can not use expressions ( ]2^J2^| ) to calculate spin distribution in the 



/-skyrmions with large effective radius. But this problem may be 
solved by taking bare radius R in the /-skyrmion equal to corre- 
sponding effective radius r\ = However, in this case we can- 
not make use of advantages provided by a simple Gaussian form 
for the overlap integral of single skyrmionic states available only 
for R w < 2 A (see (|7|)). The orbital motion with large bare radius 
acquires some unconventional properties compared to the case of 
" small orbits" . Spin distribution in these states can be well approx- 
imated by the following expressions 

2A 2 

S z (z) = -1+ (27) 
l{\ 2 + zz + R 2 ) 2 -AzzR 2 

„ . » A , zz — \ 2 — R 2 . 

S + (z) = - 1 + 28 

z ^(A 2 + zz + R 2 ) 2 - AzzR 2 

This expressions can be easily obtained if one neglects overlap of 
skyrmions in different points at all. Here, R equals to ri = as 
was explained before. Numerical calculations show that if we take 
account of the skyrmion overlap the effective radius of the state 
with bare radius R may be changed but orbital momentum can not 
strongly deviate from 27r ^ 2 fi2 because otherwise discretness effects 
become again significant. Expression (|27j) shows that maximal spin 
contraction occurs not in the origin as is the case for small skyrmions 
but close to the points with | z \— R. 
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Strictly speaking, the correct description of large skyrmion im- 
plies that one has keep both bare radius R and orbital momentum 
I to characterize it because matrix element of the Hamiltonian and 
overlap integral are nontrivial when i? 12 > A. But for large skyrmion 
radia such a description is redundant at least when we are interested 
in the properties of the zero Landau level only The corresponding 
effects would be negligible. 

As to the energy of these states its difference from the classical 
skyrmion energy is negligibe for large skyrmion radia so both large 
and small skyrmions may be said to constitute zero Landau level of 
skyrmion motion. In principle, it is possible to obtain higher Lan- 
dau levels but we think it requres detailed account of discretness 
effects and is beyond the scope of this paper. The calculation of the 
fine structure of the skyrmion energy due to translational effects 
is prevented by discretness effects in the framework of the present 
continual approximation since it does not allow to have arbitrary 
values of bare radius R for a given value of I so one can not real- 
ize variational procedure (0,0) similarly to the case of a vortex in 
superfluid 0. 

Both types of orbital skyrmions represent very interesting objects 
for further investigations. 

4 Conclusions 

Making use of the coherent state approximation allows to obtain 
an analytical form for an effective hopping integral for single stiff 
Belavin-Polyakov skyrmion in 2D isotropic ferromagnet. It appears, 
the skyrmion hopping is restricted by distances larger than its effec- 
tive diameter, thus resulting in a very unusual quasiparticle prop- 
erties. We introduce two types of orbital skyrmionic modes, the 
small and large /-skyrmions, respectively. For the former it is found 
both the spatial distribution of a single skyrmionic density, and 
the mean spin density distribution. The small /-skyrmions with 
nonzero orbital momentum are characterized by a strong quantum 
spin contraction effect. Translational symmetry of the Heisenberg 
Hamiltonian and topological Berry phase in overlap integral implies 
that there exists highly degenerate set of skyrmion states in analogy 
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with a zero Landau level. Based on the calculated expression for the 
transfer integral of the skyrmion we confirm existence of this highly 
degenerate set . 

It should be emphasized that the energy of the above addressed 
quantum small /-skyrmions with bare radius R < A is always equal 
to the classical single skyrmion energy independently both of R and 
I value. However, by forming the appropriate combinations of single 
skyrmion states centered at different space points with separation 
R12 > 2A the energy of the classical state can be lowered. Ad- 
dressing for example the combinations of two different states with 
separation between them exceeding 2A we obtain usual bonding and 
antibonding states. Energy of such a bonding state is 

AS 2 Jtt 
1- I S 12 (R) | 

and has a minimum at R12 = 2A. This implies that skyrmion states 
prefer to be separated a distance 2A from each other to lower the 
energy. Naturally, this effect will lower the energy of the orbital 
skyrmions with rather large bare radius R > A as compared to the 
energy of classical skyrmion, but for large skyrmion radia A this 
lowering is obviously negligible. Thus we may say that both large 
and small /-skyrmions form unified zero Landau level. 

The interaction of orbital skyrmions having the same orbital mo- 
mentum but different bare radius will result in a complex skyrmionic 
energy spectrum, which analysis is a challenge for a further study. 
However, continual approximation used here is inadequate for cal- 
culation of that spectrum because real lattice on which skyrmion 
moves has no true rotational symmetry. 

A further elaboration of the quasiparticle approach to skyrmions 
implies a detailed description of the internal structure of the orbital 
skyrmions without any constraint on its bare radius, and realization 
of variational procedure 0. An appropriate work is in progress. 
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